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Abstract 

fy^ In this article we use Landau- Kolmogorov inequalities to obtain new 

p I ■ inequalities between some common probability metrics. We illustrate 

the use of these inequalities by providing rates of convergencce in novel 
local limit theorems for the magnetization in the Curie- Weiss model 
at high temperature, the number of triangles and isolated vertices in 
Erdos-Renyi random graphs, as well as the independence number in a 
geometric random graph. We also obtain new rates of convergence in 
. other probability metrics for some of these examples. 
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The classical Landau inequality (see lHardy. Landau, and Littlewoodl (11935 



Section 3)) relates the norm of a function with that of its first and second 
derivatives. There are many ex tensions and embel li shme nts of this inequality 



in the analysis literature; see iKwong and Zettll (|1992l ) for a book length 
treatment. The following theorem is a discrete version of such inequalities. 
For a function / with domain the integers, denote for 1 ^ p < oo, 



1/p 



and H/ 1 |oo = sup igs and also define the operators A n recursively by 

A°f(k) = f(k) and A n+1 f(k) = A n f(k + 1) - A n f(k). 
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Theorem 1.1 (|Kwong and Zettl (| 1991 Theorem 4.1)). Let k and n be 



integers with 1 ^ k < n and let 1 ^ p, q, r ^ oo given. There is a positive 
number C such that 

l|A fc /II^C H/HJ-" \\A n ff r 
for all f : % — > IR, with \\f\\ p < oo and \\A n f\\ r < oo, where 

k — 1/ q + 1/p 



n — 1/r + 1/p' 



if and only if 



n n — k k 

- s$ + -. 

q p r 



Remark 1.1. Much of the literature surrounding these so-called Landau- 
Kolmogorov inequalities is concerned with finding the optimal value of the 
constant C. In the case that n = 2 and either p = q = r = 1, or n = 3 , 
p = q = oo, and r = 1, we can take C = y/2; see iKwong and Zettl (jl992l . 



Theorem 4.2). These are two of the main cases discussed below. Also, an 
inductive argument in n implies that in the former case above we may take 
C = 2( n " 1 )/ 2 for n ^ 2 and in the latter C = 2( n " 2 V 2 for n ^ 3. These 
facts are not critical in what follows, so for the sake of simplicity we do not 
discuss such constants in further detail. 

Applying Theorem 11.11 to the case / = F — G and k = r = 1, and using 
the triangle inequality we obtain the following result. 

Corollary 1.2. Let n > 1 be an integer and 1 ^ p, q ^ oo be such that 

n^n-l+p^ 
q p 

If F and G are functions with domain the integers, then there is a positive 
number C depending only on n, p and q such that 

|| AF - AG\\ q ^C\\F- G\\l~ (|| A n F||i + HA^Hi/, 

where 

P= 1 - l ' q + l l P . (1.2) 
H n-l + l/p y ' 
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The importance of Corollary 11,21 stems from the fact that if F and G are 
distribution functions of integer supported probability distributions, then 
some well known probability metrics can be expressed as 



d K (F,G) Hl^-Glloo 
d w (F,G) HI* 1 -GHi 
d loc (F,G) =||AF-AG|| K 
d TY (F,G) =±\\AF- AG\\ 
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(Kolmogorov metric), 
(Wasserstein metric), 
(local metric), 
(total variation metric). 



Note that d loc obtains its name from the fact that it is intimately connected 
to local limit theorems; see Section [H 

Corollary 11.21 implies, for example, that 



In this article we think of F as being a distribution of interest (e.g. the 
number of triangles in a random graph model) and of G as a well-known 
distribution which we are using to approximate F (e.g. a discretized normal 
or a translated Poisson distribution). Thus, if a bound on djz(F, G) is known, 
then in order to bound dy oc (F,G), it is enough to obtain good bounds on 
||A 3 F||i and ||A 3 G||i. Whereas it will be easy to bound ||A 3 G||i, since G is 
assumed to be an explicit distribution, estimating ||A 3 F||i is more involved. 
Therefore, (I1.3P is only of use if we have techniques to bound ||A 3 F||i in 
non-trivial situations. 

The purpose of this paper is two-fold. Firstly, in Section [2] we explore 
inequalities of the form (|1 .31) : this is more or less a straightforward conse- 
quence of Corollary 11.21 Secondly, in Section [3] we develop tools to bound 
||A 3 i ? ||i and the analogous quantities appearing on the right hand side of 
generalizations of (|1.3p . In Section UJ we illustrate our approach in a few 
applications, in particular we obtain new local limit theorems with rates 
for the independence number of a geometric random graph, for the num- 
ber of isolated vertices and triangles in Erdos-Renyi random graphs and the 
magnetization in the Curie- Weiss model. We also obtain other new limit 
theorems and rates for some of these applications. 



The next result is our main theorem, and is a slight generalisation of Corol- 
lary 11.21 The main purpose of the generalisation is to handle certain av- 
eraged local limit theorems where the corresponding pointwise results may 
not hold. 



d loc (F,G)^Cd K (F,G) 



^(HA^IIi + IIA^IIi) 



1/2 



(1.3) 



2 MAIN RESULT 
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Let us introduce a last bit of notation before stating the main result. Let 
F be a distribution function with support on 7L. If m is a positive integer, 
define 

pm{j) = FU) + --- + F(j-m+l) _ 
m 

this is the distribution function of the convolution of F with the uniform 
distribution on {0, . . . , m — 1}; note that F l = F. 

Theorem 2.1. If I 1 and rn 1 are integers, then there is a constant 
C > such that, for all distribution functions F and G of integer supported 
probability distributions, 

di{F m ,G m ) < Cd 2 (i ? ,G) 1 ~ /3 (||A m F m ||i + ||A m G m ||i) /3 
for the following combinations of d\, di and f3: 





di 


d 2 


P 


(0 


^loc 


C?TV 


l/l 


(it) 


d loc 


d K 


l/l 


{Hi) 


^loc 


c?w 


2/(1 + 1) 


(iv) 




c?w 


i/G + i) 


(v) 


r/ K 


dw 


i/a + i) 



Proof. To prove (ii)-(iv), apply Corollary 11.21 to the functions F m and G m , 
n = I + 1, and with the following values of p and (/: 

(h) q = oo,p = oo, (Hi) q = oo,p=l, (iv) q = l,p=l; 

then use d2(F m ,G m ) ^ d2(F,G). For (i) and (iu) use (ii) and (u), respec- 
tively, and then use the fact that c?k ^ ^tv- D 

Let us make a couple of remarks at this point. First, we mainly use 
Theorem 12.11 with m = 1, where its meaning is most transparent. For 
m > 1, the following easy lemma should clarify our use of Theorem 12.11 
below. 

Lemma 2.2. If X andY are integer valued random variables with respective 
distribution functions F and G, then 

sup|P[X G (k,k + m}] -P[y G (jfe,Jfe + m]]| = md loc (F m ,G m ). 
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Second, Theorem 12.11 is really a special case of Theorem 11.11 with k = 
r = 1, but it is clear that similar statements hold by applying Theorem ll.il 
to other values of k and r. We choose the value r = 1 because we are able 
to bound ||A n F||i. Using the obvious inequality ||A n F|| 00 < ||A n+1 F||i, 
we could also usefully apply Theorem 1 1 . 1 1 wit h r = oo, but this change has 
no effect on the value of [3 for a given k, q and p. The term HA 2 ^^ also 
appears in the local limit theorem results of McDonald! ( 1979) and Davis and 



McDonaldll995|). However, the crucial advantage of ||A n F||i over ||A n F| 



is that the former is — as we will show — amenable to bounds via probabilistic 
techniques, whereas the latter seems difficult to handle directly. 

3 ESTIMATING THE MEASURE OF SMOOTHNESS 

In this section, we develop techniques to bound ||A n F m ||i. To lighten the 
notation somewhat, write 

D n , m (F) :=m||A n+1 F m ||i, 

and D n for D n ^\. Furthermore, define recursively the difference operators 

A n m F{j) = A^F(j + m) - A^F(j), 

where A° m F(j) = F(j). 

The key to bounding D n ^ m (F) is the following representation. 

Lemma 3.1. Let n and m be non-negative integers. Let W be a random 
variable with integer support. Then 

D n>m (J?(W))= sup EA^g(W). 

Proof. We only show the case n = 1; the general case is similar. Denote 
Pi = P[W = i] for i £ 7L. We have 



EA m9 (U/) = ^2 Pi {g(i + m) - gii)) 

= ^2(Pi-rn ~ Pi)g{i) 

E( Pi-ra H h Pi-i Pi 
— m+1 
sW 

^{^F m (i-2)-AF m (i-l))g{i) 
^A 2 F m ( l -2)g( l ). 



m 

■iGS 



-m 
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We see that for all g such that WgW^ sC 1, E,A m g(W) m||A 2 F m ||i, and 
choosing g{i) = — sgn A 2 F m (i — 2) implies the claim. □ 

The following sequence of lemmas provide tools for bounding D n ^ m {F). 
The proofs are mostly straightforward. We assume that all random variables 
are integer valued. 

Lemma 3.2. Let n and m be positive integers. Let W be a random variable. 
Then 

D n , m {&(W)) ^ mD nA (J?(W)). 

Proof. Using a telescoping sum argument, this is immediate from the defi- 
nition. □ 

Lemma 3.3. Let n and m be positive integers. Let W be a random variable, 
and let T be a a -algebra. Then, 

D n , m (Sf(W)) ^ ®D n>m {J?(W\T)). 

Proof. If / is a bounded function, then 

|EA£/(W)| < miKJ(W)\T)]\ < ll/IUED^J^WJ-)). □ 

Lemma 3.4. If Xi and X 2 are independent random variables, then, for all 
ni,n 2 ,m > 1, 

D ni+n2 , m {^{X l + X 2 )) D num (^'(X 1 ))D n2 , m (^'(X 2 )). (3.1) 

If Xi, . . . ,Xjy is a sequence of independent random variables and n ^ N, 

n 

£>„, m (JSf (X 1 + ■■■ + X N )) < [] D lirn {#{Xi)). (3.2) 

i=l 

Proof. Let / be a bounded function. Define g{x) := EAJ^ 2 /(s + X 2 ) and 
note that, because of independence, 

EA£+*V(Xi + X 2 ) = EA^pd). 

Hence 

(jZ?(Xl +X 2 )) < -Dni,m(^l)lblloo < Aii,m(^l)Ai 2 ,m(^2)||/||oo 

which proves (|3.ip ; (|3.2|) follows by induction and Lemma 13.31 □ 



1. 



Lemma 3.5 (jMattner and Rood (|2007l ). Corollary 1.6). Let X u X 2 ,...,X N 
be a sequence of integer-valued random variables and Sn = Si=i -^i- Then 

i=l 

The following two theorems are our main contributions in this section. 

Theorem 3.6. Let (X,X') be an exchangeable pair and let W := W(X) 
and W' := W(X') take values on the integers. Define 

Qm(x) = P[W' = W + m\X = x] 

and q m = EQ m (X) = P[W' = W + m]. Then, for every positive integer m, 

Proof. To prove the first assertion, we must bound |EA m <?(W)| for all g 
with norm no greater than one. To this end, exchangeability implies that 
for all bounded functions g 

= q^{l[W = W + m]g(W) - l[W = W - m]g(W)} 
= q- l ^{Q m {X)g{W + m) - Q_ m (X)g(W)}, 

so that 



\^A m g(W)\ 

= |E{(1 - q^Q m {X)g{W + k) - (1 - q- 1 Q. m {X)g{W)}\ 



jM{g(W + m) 2 } Var Q m (X) + ^{g{W) 2 } Var Q- m (X) 

q m 



(3.4) 



where in the inequality we use first the triangle inequality and then Cauchy- 
Schwarz. Taking the supremum over g with ||g||oo ^ 1 in (|3.4p proves the 
theorem. □ 

Theorem 3.7. Let (X, X',X") be three consecutive steps of a reversible 
Markov chain in equilibrium. Let W and W' be as in Theorem \3.(A and, in 
addition, W" := W{X"). Define 

Qm um2 (x) = P[W' = W + mi, W" = W + m 2 \X = x]. 
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Then, for every positive integer m, 

D 2 , m {J?(W)) < 4r (2 Var Q m (X) + E|Q m>m (X) - Q m (X) 2 \ 
+ 2VarQ_ m (X) + E|Q 

— m,— m 

Proof. Similar to the proof of Theorem 13.61 we want to bound m/S? m g{W) 
for all 5 with norm no greater than one. We begin with the trivial equality 

E{I[W = W + m, W" = W + m]g(W + m)} (3.5) 
= E{I[W' = W + m,W" = W' + m]g{W')}. (3.6) 

Conditioning on X in ()3.5[) and on X' in (|3.6p . the Markov property and 
reversibility imply 

^{Q m ,m{X)g(W + m)}= M{Q m (X)Q_ m (X)g(W)}, 

and similarly 

E{Q 

—m,—m 

Using these two equalities coupled with (|3.3[) it is not hard to see that for 
bounded g 

E,A 2 g(W) 

= M{g(W + 2)((1 - q^Q^X)) 2 + g- 2 (g m , m (X) - Q m {X) 2 )) } 
- 2E<[ 5 (^ + 1) (1 - q^QmiX)) (1 - g- 1 Q_ m (X)) } 
+ e{ 5 (W)((1 - q-iQ-niX)) 2 + q- 2 {Q- m ,- m (X) - Q- m (X) 2 )) }. 

The theorem now follows by taking the supremum over g with norm no 
greater than one and applying the triangle inequality and Cauchy-Schwarz. 

□ 

4 APPLICATIONS 

Because we are going to work in the total variation and local limit met- 
rics, we need to use a discrete analog of the normal distribution. For the 
sake of concreteness, we use the translated Poisson distribution. We say 
that the random variable Z has the translated Poisson distribution, denoted 
Z ~ TP(/i, a 2 ), if Z - [p. - a 2 \ ~ Po(<r 2 + 7), where 7 = y. - a 2 - [p, - a 2 \ . 
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Note that 1EZ = fi and a 2 ^ VarZ ^ a 2 + 1. The following lemma es- 
sentially states that we can use the translated Poisson distribution as a 
discrete substitute for the normal distribution and also provides bounds on 
the appropriate smoothing term. 

Lemma 4.1. I//i6E and a 2 > 0, then as a — > oo, 

D k , m (TF(pL,a 2 )) = 0(0, (4.1) 
d K (TP(/i,a 2 ),N(^,a 2 )) =0(0, (4-2) 
d w (TP(/0,N(/0) =0(1), (4.3) 



and 



sup 

fees 



TP(/i,a){fc}-^==exp 



0(0- (4-4) 



Proof. Since Po(a 2 + 7 ) can be represented as the convolution of [o~ 2 \ Poisson 
distributions with means at least 1, (14. ip follows from Lemmas 13.21 13.41 
and 13.51 Using the same decomposition, the remaining properties follow 
fr om the classica l theory of sums of independent random variables; we refer 



to lPetrovl (|l975h . □ 



Before moving on to our applications, we make a few remarks about 
Lemma [4. H and its use in what follows. First note that the rates obtained in 
Lemma [4. II hold in general for sums X\ + • • • + X n of independent identically 
distributed random variables with integer support and Di(Xi) < 2. Also, in 
order to appreciate the Wasserstein bound in Lemma 14. 11 the reader should 
keep in mind that for random variables X and Y and any positive constant c, 
we have 

d w (^{cX),^{cY)) =cd w (^{X),^{Y)). 

Finally, equation (|4.4p is just the local limit theorem for the translated 
Poisson distribution. Such a statement is only informative if the right hand 
side of (j4.4p is o(<r~ 1 ), because the left hand side is trivially O^a" 1 ). In this 
section we will prove bounds for d loc ( = 5f(VF),TP( / u,c 2 )), which are better 
than 0(a~ 1 ) and therefore, by means of (|4.4p . will lead to a local limit 
theorem for W along with a rate of convergence. 



4.1 Embedded sum of independent random variables 

As the first application of our general theory, we consider the case where W 
has an embedded sum of independent random variables. This setting has 
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been the most prominent way to prove local limit th e orems by prob a bilisti c 
argumen t s; see , for example, bavis and McDonaldl (119951 s ). IRollinI (j2005h . 
Barbour! d2009h. IBehrisch. Coi a-Oghl an. and Kand and Penrose anc 

Peres (|2011 ). In this case, our theory can be used to obtain rates for the 
accuracy of the approximation if rates of convergence for W to its limit are 
known. 

Let W is an integer-valued random variable with variance a 2 and let T 
be some c-algebra. Assume that W allows for a decomposition of the form 



N 



W = Y + Y,Zi 



(4.5) 



i=l 



where N is ^-measurable, and where, conditional on J 7 , we have that 
Y, Z\, . . . , Zn are all independent of each other. Note that in what follows, 
the distribution of Y is not relevant. 

Theorem 4.2. Let W = W(o~) be a family of integer-valued random variable 
with Variy = a 2 satisfying (|4.5p . Assume there are constants u and (3, 
independent of a 2 , such that, conditional on J- , 



< u < 1 - \D x (Z t ) 
for all 1 N, and such that 

P[N < Pa 2 } = 0(d- fc ) 
as a — > 00 for some k ^ 2. With W = (W — EW)/<r, then as a —> 00, 



(4.6) 



(4.7) 



d loc {j?(W),TP(VW,a 2 )) =01 



a 



and 



d loc (j?(W),TP(EW,a 2 )) =0 



l-2/(fc+l)- 



a 



Retaining the previous hypotheses, if (|4.7p ZioWs now /or some ^ 1, £/ien 
d Ty (^(W),TP{]EW,a 2 )) = 0(d w (^(W),<Z>) k/(k+iy 



as a — > 00. 
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Proof. First, consider the setup conditional on T . Divide the sum Z\ + • • • + 
Zjy into k successive blocks, each of size [N/k\ , plus one last block with less 
than [N/k\ elements, which we will ignore in the following estimates. By 
Lemma 13.41 we have 



D k {W\F) ^f[m, 



1=1 



where 



Vi = Dx{ z {i-\)\_N/k\+i H + z i[N/k\) 

/ l\N/k\ v-1/2 

\ i=(l~l)[N/k]+l J 

where we have used Lemma 1331 From (|4.6p . for I = 1, . . . , k, 

m < 2 A 



7rLiV/A;Jny ' 

Now, assume without loss of generality that a 2 > k//3. Since we can always 
bound Dk(W) by 2 fc , we have 

D k (W\T) ^ 2 fc I[iV < P<j 2 ) + 21[N > Pa 2 } 



TTU{(3CT 2 — k) 

Hence, by Lemma 13.31 

8k ^ k/2 



D k {W) < MD k {W\F) ^ 2 k ¥[N < Pa 2 } + 2 



iru(f3a 2 — k) 



which is 0(a~ k ) as a — > oo. After noting that since W is integer valued, the 
scaling a" 1 = 0(dy/(j!?(W), $) A {J£{W\ ) , the claims now follows 
easily from (ii), (in), and (iv) of Theorem 12.11 and Lemma 14.11 □ 

Note that, under the stated conditions, Theorem 14.21 implies the LLT 
for W if it satisfies the CLT, as the latter also implies convergence in the 
Kolmogorov metric. If a rate of convergence is available, Theorem 14.21 also 
yields a rate of convergence for the LLT. 

Let us apply now Theorem 14.21 to the so-called independence number of 
a random graph. The independence number of a graph G is defined to be 
the maximal number of vertices that one can pick from the graph so that 
no two of these vertices are connected. 
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Consider the foll owing random graph model , which is a simplified version 



of one discussed bv lPenrose and Yukichl (|2005l ). Consider an open set U C 
R d , d ^ 1, of finite volume, which, without loss of generality, we assume to 
be 1. Let X be a homogeneous Poisson point process on U with intensity 
A with respect to the Lebesgue measure. Define G(X,b) to be the graph 
on the vertex set X by connecting two vertices whenever they are at most 
distant r apart from each other. In the context of this random geometric 
graph, the independence number is the maximal number of closed balls of 
radius r/2 with centers chosen from X, so that no two balls are intersecting. 

Theorem 4.3 ([Penrose and Yukichl ([20051 1). Let W b = /3(G(X, bX' 1 ^)) . 



Then, for b small enough, we have Var W b X A and 
as A — » oo. 

Note that the condition "for b small enoug h" is necessary to guarantee 



the asymptotic order of the variance of W b ; see Penrose and Yukichl (]2005l ) 



for more details. We can give a local limit result as follows. 



Theorem 4.4. Under the assumptions of Theorem \4-3j we have that for 
every e > 0, 

d loc (j?(W b ), TP (TEW b ,VarW b )) = 0(A^ 1+£ ), 

as A — >■ oo 

Proof. Let r = b\~ 1 / d /2. Denote by B r (x) the closed ball with radius r and 
center x; define dB r (x) = B2 r (x) \ B r (x). Now, choose n non-intersecting 
balls in U, each of radius 3r and centers x\,...,x n ; it is clear that it is 
possible to have n x A. For ball B^ r (xi), define the indicators 

Ii = l[dB r (xi) n X is empty], Ji = l[B r (xi) n X is not empty]. 

Note that the Ii are independent and identically distributed and, hence, 
N = Yji=i^i ~ Bi(n,p) with p = E7j being bounded away from 1 and 
as A — > oo. We let T = <j(Ji, . . . ,I n )- Furthermore note that if I\ = 1, 
then Ji is exactly the contribution of the ball B2 r (xi) to the independence 
number W b , as within B r {xi) all the vertices are connected and there is no 
connection to any other vertices outside B2 r {xi). Therefore we can find Y 
such that 

N 



W n = Y + Y,JK j , 



12 



where K\, . . . ,Kjy are the indices of those balls with Ij = 1. Given 
note that Jk 1 , ■ ■ ■ , Jk n are independent Be(g), with with q = EJ^ being 
bounded away from and 1, and they are also independent of Y . This 
implies condition (|4.6|) for u = q A (1 — q) which is bounded away from as 
A — > oo. Using usual exponential tail bounds for the binomial distribution, 
it is easy to see that, for every k, one can find /3 such that (|4.7p holds. In 
combination with Theorem 14.31 Theorem 14.21 now yields the claim. □ 

4.2 Magnetization in the Curie- Weiss Model 

Let f3 > 0, h S H and for a G { — 1, l} n define the Gibbs measure 



P(<t) = Z- 1 exp|-^fj J cj j + /i^ail 



(4.8) 



where Z is the appropriate normalizing constant. 

This probability model is referred to as the Curie- Weiss model and a 
quantity o f inter est is the magnetization W = ^ • cr, of the system. See 



Talagrandl (|201ll ) for an introduction to the vast literature around these 



models. We have the following known results. 



Theorem 4.5 ( Ellis. Newman, and Rosen ( 1980l . Theorem 2.2)). Let a have 



law given by (|4.8|) with < f3 < 1, and fteE. ^Uso, Zei mo &e the unique 
solution to 

m = tanh(/3m + h). 
Then there is a constant C depending only on (3 and h such that, as n — >■ oo ; 

d K (^(^),N(0, I ^,))^0. 



Theorem 4.6 (jEichelsbacher and Lowd (j2010l . Theorems 3.3 and 3.7)). Let 



a have law given by (|4.8p wii/i < /3 < 1, anrf h = 0. Then there is a 
constant C depending only on f3 such that 

d K {^(n-^ 2 W),N(0,(l - P)- 1 )) ^Cn- 1 ' 2 ; 

the same bound holds for the Wasserstein metric. 

Let us now give a total variation metric and local limit version of these 
results. 
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Theorem 4.7. Let W be as before and let 5 = 6(n) = (1 — (— l) n )/2. For 
< f3 < 1 and h = 0, there is a constant C that only depends on (3 such 
that 

d }oc (3?((W + 5)/2),TP(0, 3p ^ y )) < Cn^'\ 
d TY {j?((W + <y)/2),TP(0, ^)) < Cn-V3. 

If < (3 < 1, fteE, and mo is as in Theorem \4.5[ then there is a constant 
C that only depends on f3 and h such that as n — > oo, 

A ( q> ( W+5 \ Tp /m n(l-mg) \\ _ / -l/2\ 

Proof. The proof follows easily from (it), (Hi) and (iv) of Theorem 12.11 with 
m = 1 and Z = 2, Lemma 14.11 Theorems 14.61 and 14.51 and the bounds on the 
smoothing terms read from Lemma 14.81 below. □ 

Remark 4.1. In the critical case where /3 = 1, a bound is known on the 
Wasserstein distance between the magn etization (appropriately norma lized) 
and its non-normal limiting distribution Eichelsbacher and Lowej ( 201(1 The- 



orem 3.3). The smoothing bounds of Lemma 14.81 below apply to this case 
with h = (see Remark 14. 2p and could be used to prove new bounds on 
the total variation distance between the magnetization and a discrete ver- 
sion of this limiting distribution. However, we omit this result due to the 
inappropriate amount of space it would take for a precise formulation. 

In order to apply our theory we use Theorem 13. 71 to bound the necessary 
smoothing terms. For this purpose, let a as in the theorem and a' be a 
step from a in the following reversible Markov chain: at each step of the 
chain a site from the n possible sites is chosen uniformly at random and 
then the spin at that site is resampled according to the Gibbs measure (|4.8j) 
conditional on the value of the spins at all other sites. Let W = Y17=i °» 
and W = J27=i °'i an d note that (W, W) is an exchangeable pair. Finally, 
define 

Q m = P[W' = W + m\a], 

q m = EiQ m , and 

Q mi , m2 = P[W" = W + mi, W" = W' + m 2 \a}. 

where W" is obtained from W' in the same way that W is obtained from 
W (i.e. (W,W' ,W") are the magnetizations in three consecutive steps in 
the stationary Markov chain described above). We will show the following 
result. 
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Lemma 4.8. Let < ft < 1, /iGR and m = — Y2™=i a i? an d let m o be the 

unique solution to 

m = tanh(/3m + h). 



Then, for k = ±2, 



1 — m, 



^ C ( \m — mo | + 



and 



1 — m 



Qk 



\Qk,k ~ Qt\ = CKn" 1 ), 



0(n- 1 ' 2 ), Var(Q fc ) = 0(n~^ 



(4.9) 

(4.10) 
(4.11) 

(4.12) 



Proof. We only consider k = 2, the case /c = —2 being similar. An easy 
calculation shows that, 



P(^ = l|(^) jV 



so that for m 8 := n~ l Ylj^i a ii we have 

1 1 — <7j tanh(/3m,j + /t) + 1 



<3s 



n ^ 2 



Now, some simplification shows 

1 m tanh(/3m + /i) 
Q2 = 4"T + 4 ' 



1 — m) 



1 

H y^(l - o-j)(tanh(/3mj + h) - tanh(/3m + h)). 



i=l 



Thus we find 



Qi 



1 — m, 



^ -|m — mo| + -|tanh(/3mo + /t) — tanh(/3m + /i)| 

+ — |mtanh(/3m + /i) — motanh(/3mo + /i)| 
1 ™ 

+ — VYl - o-j)| tanh(/3mj + /i) - tanh(/3m + 



i=l 
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Since tanh(x) G (—1,1) is 1-Lipschitz and — 1 ^ m ^ 1 the first part of the 
claim now easily follows. 

For the second assertion, note that 



Q2,2 = -^2 Y,^ 1 ~ ^)(tanh(/3mi + h) + l) 



x (1 - cr J )(tanh( / 5m i) j + /3n _1 + h) + 1), 
where rriij = /3n _1 Ylk^i j a ki an d also that 

Q2 = ttt - 2 X^ 1 ~ o-OCtanhC 9 "*! + M + - ^)(tanh(/3?7ij + /i) + 1). 



1,3 



We can now find 

\Q%1 - Oil < r-2 X^ 1 " ^)(tanh(/3mi + h) + 1) 



+ ^ Z^ 1 ~ ^X 1 ~ tanh(/3mi + h) + 1| 



x |tanh(/3mjj + /3n 1 + h) — tanh(/3mj + 



Straightforward estimates now yield (|4,lip . 

The remaining assertions (|4.12p follow from (|4.10p and E|m — mo\ k = 
0(n~ k / 2 ) which is obtained from standard concentration results; see e.g. 
Chatter i eel (|2007l . Proposition 1.3). □ 



Remark 4.2. Lemma 14.81 easily extends to the case (3 ^ 1 if h = as the 
solutions to (14.9P are symmetric around and therefore ttt-q is the same for 
any such solution. 

4.3 Isolated vertices in the Erdos-Renyi random graph 

In this and the next section we will derive LLTs for the number of isolated 
vertices and triangles in the Erdos-Renyi random graph, respectively. LLTs 
related to random gr aphs using probabilistic methods seem to be sparse. 
Behrisch et~ai1 (j2O10h obtain LLTs for the size of the maximal component 
in hypergraphs using an embedding idea related to Secti on [4TT1 Similar re- 
sults on the maximal component such as lStepanov (Il970h are more analytic 
i n nat ure. Using combinatorial arguments, Bender. Canfield. and McKay 
(|l997l ) obtain very detailed information about the number of isolated ver- 
tices for graphs with a given number of vertices and edges. To the best of 
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our knowledge, the following results on the number of isolated vertices and 
number of triangles are new. 

In order to estimate the smoothness in the random graph examples we 
will make use of Theorems 13.61 and 13.71 and then combine them with Theo- 
rem [2D] to obtain the local limit results. Define G = G(n,p) to be a random 
graph with n vertices where each edge appears with probability p, indepen- 
dent of all other edges. From this point, we have the following theorem. 

Theorem 4.9 ( Kordecki ( 199Cj| ); Barbour. Karohski. and Rucihski ( 19891 )). 

Let W = W(n,p) be the number of isolated vertices of G(n,p), and let W 
be W normalized to have zero mean and unit variance. Then W converges 
in distribution to the standard normal if and only if 



lim n p = lim (log(n) — np) = oo. (4-13) 

n— >-oo n— >oo 

In that case, with a\ = Var W , 

d K {J?{W)^) = 0(o- 1 ). 



The conditions of convergence was proved by iBarbour et al.l (11 989). 
where as the bounds for the Kolmogorov metric was obtained by iKordecki 
(jl99(]l ). 

We now summarize the results derived from the combination of Theo- 
rem 14.91 and Theorem 12.11 For two distribution functions F and G with 
integer support, let 

d&(F,G) = ||AF m -AG m ||oo. 

Note that d loc = and recall also the inequality given by Lemma 12.21 

Theorem 4.10. Let W = W(n,p) be the number of isolated vertices in an 
Erdos-Renyi random graph G(n,p) and W be W normalized to have zero 
mean and unit variance. With fi n = WW and a\ = Var(M^), we have the 
following. 

(i) If lim (log(n)— np) = oo, and either lim np = oo or lim np = c > 0, 

n— >oo n— >oo n— >oo 

then <j2 x ne~ np and 

d loc (j?(W),TP(n n ,a 2 n )) =0(<7- 3 / 2 ). 
(ii) If lim np = and lim n 2 p = oo, then a 2 x n 2 p and 

n— ¥oo n— >oo 

d loc (j?(W),TF(pL n ,al)) = Ofc-^np 3 / 4 )- 1 ), 
d 2 oc {J?(W),TF(n n ,c- 2 n )) =0{a^ 2 ). 
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Proof. After noting 



a 2 n = n{\ - p) n ~ l [1 + (np - 1)(1 - p) n ~ 2 ] , 

the result follows from (n) of Theorem 12.11 with 1 = 2, using the known rates 
stated above in Theorem 14.91 coupled with Lemma I4.1| and bounds on the 
smoothing quantities provided by Lemma 14.111 below. □ 

Remark 4.3. The bounds on the smoothing quantities presented below can 
be written in terms of n and p, so that the asymptotic results of the theorem 
can be written explicitly whenever the bounds on d^(W, <3?) are also explicit. 

Remark 4.4. The second case of the theorem is interesting and deserves 
elaboration. In some regimes, our bounds do not imply a LLT in the natural 
lattice of span one (e.g. p x n~ a , where a > 4/3), but we can obtain a useful 
rate in the d 2 oc distance. This implies that the approximation is better 
by averaging the probability mass function of W over two neighbouring 
integers and then comparing this value to its analog for the normal density. 
One explanation for this phenomenon is that in such a regime, the graph 
G(n,p) will be extremely sparse so that parity of W will be dominated by 
the number of isolated edges. In other words, up to a coin toss with small 
success probability, n — W will be approximately equal to twice the number 
of isolated edges, in which case we would not expect the normal density to 
be a good approximation for each point on the integer lattice. 

Lemma 4.11. Let W = W(n,p) be the number of isolated vertices in an 
Erdos-Renyi graph G(n,p). 

(i) If lim (log(n) — np) = oo, and either lim np = oo or lim np = c > 0, 

n—too n— >oo n— >oo 

then 

D 1 (W) = 0(a- 1 ), D 2 (W) = 0{o-~ 2 ). 

(ii) If lim ra ^. 0O np = and linin^oo n 2 p = oo, then 

D 1 (W) = 0((^„)" 1 ) ) D 2 (W) = 0((np)- 1 (7- 2 ), 
D li2 (W) = O^ 1 ), D 2)2 (W) = 0{a~ 2 ). 

Proof. We will apply Theorems 13.61 and 13.71 by constructing a Markov chain 
on graphs with n vertices which is reversible with respect to the law of 
G(n,p). From a given graph G, define a step in the chain to G' by choosing 
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two vertices of G uniformly at random and independently resampling the 
"edge" between them. It is clear that this Markov chain is reversible with 
respect to the distribution of G(n,p). 

Now, let W = W(n,p) be the number of isolated vertices of the Erdos- 
Renyi graph G = G(n,p) and and W be the number of isolated vertices 
of G' , a step in the chain described above. In order to compute the terms 
needed to apply Theorems 13.61 and I3.7( we need some auxiliary random 
variables. Let Wk be the number of vertices of degree k in G (so that 
Wo = W), and E 2 be the number of connected pairs of vertices each having 
degree one (that is, E 2 is the number of isolated edges). In the notation of 
Theorems 13.61 and 13.7} we have 

(W 1 - 2E 2 ) W(n - W) 
Qi(G) = (1-P), Q-i(G) = ^ p, 

9 (Wi-2E2\ A(W\(n-W+l\ 

QiAG) = 1 2 U i-P) 2 , Q-xMG) = Uj r L 2 V , 

E ( w ) 
g 2 (G) = -^(l-p), Q_ 2 (G) = ^jp, 



^, 2 (G) = ^/(l - P)\ Q-2MG) = [ -^rf 1 P 2 - 

\2/ \2/ 



These equalities are obtained through straightforward considerations. For 
example, in order for one step in the chain to increase the number of isolated 
vertices by one, an edge of G must be chosen that has exactly one end vertex 
of degree one, and then must be removed upon resampling. For one step 
in the chain to decrease the number of isolated vertices by one, an isolated 
vertex must be connected to a vertex with positive degree. 

From this point, the lemma will follow after computing the pertinent 
moment information needed to apply Theorems 13.61 and 13.71 By considering 
appropriate indicator functions, it is an elementary combinatorial exercise 
to obtain 

EW? = 2^(1 -p)- 2 + 2 P (fj [(1 -p) 2 - 4 +p(n - 2) 2 (1 -pf n ~% 

veI = QMi-ri^ + eQp^-P) 4 "- 12 , 
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= -P) 2n " 4 [(" - 2)(n - 3)p(l -p) n - 

which will yield the results for negative jumps, and 

ew = n(i - P ) n ~\ ew 2 = n(l -p) 12 - 1 + 2( '")(l-p) 



Elf 3 = n(l - p) n ~ l + 6^) (1 " P) 2n ~ 3 + 6 ( o ) C 1 " P) 3 "~ 
EW 4 = n(l - p) 71 " 1 + 14 f ™ ) (1 - p 



+ 36(M (1-^^ + 24(^1(1-^ 

which will yield the results for the positive jumps. Theorems 13.61 and 13.71 
now give the desired rates. As an example of these calculations, note that 
for qi = 1EQi(G), we have 

Var Q X {G) _ EW? - 4EWi£ 2 + 4E#| - (EWi - 2E£ 2 ) 2 
gf ~ (EWi - 2E£ 2 ) 2 ' 

which after the dust settles is 0(n _1 e np ) in case (i) of the theorem. Simi- 
larly, since q\ = EQ_i, we have 

Var Q-i(G) _ EW 4 - 2nEW 3 + n 2 EW 2 - (nEW - JEW 2 ) 2 

q} ~ (nEW - Er) 2 ' ( ' 15) 

which is again 0(n _1 e np ) in case (i). Theorem 13.61 implies that D\(W) 
is bounded above by the sum of the square roots of the terms in (|4.14p 
and (I4.15P so that in case (i), 

D 1 (W) = 0(n- 1 / 2 e^=0(a- 1 ). 

For the second part of (i), note that 

E|Qi,i-Q?| EWi + 2ES 2 



ql ^ (EWi - 2E£ 2 ) 2 ' 
which is 0(n _2 p _1 e np ) in case (i) and 

E|Q-i,-i - Q 2 J _ E(W(n - W)|n - 2W + 1|) n + 1 



qf (nEW - E^ 2 ) 2 ^ nEW - EW 2 ' 

which is 0(n _1 e np ) in case (i), so that we have 

D 2 (W) = 0(n~ l e np ) =0(a~ 2 ). 
This proves (i); the remaining bounds are similar and hence omitted. □ 
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4.4 Triangles in the Erdos-Renyi random graph 

In this section we use Theorems 12.11 [3T6l and l3.7l to first obtain a new rate of 
convergence in the total variation distance between the normal distribution 
and the number of triangles in an Erdos-Renyi random graph. We then 
use this new rate to obtain a local limit theorem for this example. As in 
Subsection 14.31 define G = G(n,p) to be a random graph with n vertices 
where each edge appears with probability p, independent of all other edges. 
Prom this point, we have the following theorem. 



Theorem 4.12 (|Rucihskil (|l988h . iBarbour et all (|l989M . Let W = W(n,p) 



be the number of triangles of G(n,p), and let W be W normalized to have 
zero mean and unit variance. Then W converges to the standard normal if 
and only if 

lim np = lim n 2 (l — p) = oo. 

n— >oo n— ¥oo 

In that case, with a 2 = Var W, 

d w (J?(W),$) = 0(o-- 1 ). 

We now summarize the results derived from the bound of Theorem 14. 121 
coupled with our theory above. 

Theorem 4.13. Let W = W(n,p) be the number of triangles in an Erdos- 
Renyi random graph G(n,p). If n a p — > c > with 1/2 ^ a < 1 then, with 
fj, n = EW and a\ := Vav(W), we have a\ x n 3 p 3 and 

d TV (^(W),TP(fi n ,a 2 n )) = O^ 1 -")), 

and 

d loc {^(W),TP(f, n ,a 2 n )) = O^n-V-^ 2 ). 
Proof. After noting that 



a 2 = (p\l -V) + 3 (n- 3)p & (l -p)), 

the result follows from (iv) and then (i) (or {Hi)) of Theorem 12.11 with 
/ = 2 and m = 1 using the known rates stated above in Theorem 14.121 
coupled with Lemma l4.1| and bounds on the smoothing quantities provided 
by Lemma 14. 141 below. □ 
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Remark 4.5. It is worthwhile noting that we obtain the LLT only for those 
values of a for which we have JEW x Var W. In contrast, if < a < 1/2, we 
have that ~EW x n 3 ~ 3a , whereas Var TV x n 4_5 ° » EW. It is not clear if 
a standard LLT will hold in this regime as the probability mass is scattered 
over an interval that is relatively large compared to the expected number 
of triangles due to the strong dependence. The point probabilities may 
therefore behave in a much less controlled way; c.f. Remark 14.41 following 
Theorem 

Lemma 4.14. Under the assumptions of Theorem \4-13[ 

Dl (W) = 0{a~ l ), D 2 (W) = 0{a- 2 ). (4.16) 

In order to prove Lemma 14.141 we wm a PPly Theorems 13.61 and 13.71 by 
constructing a Markov chain on graphs with n vertices which is reversible 
with respect to the law of G(n,p). From a given graph G, define a step 
in the chain to G' by choosing two vertices of G uniformly at random and 
independently resampling the "edge" between them. It is clear that this 
Markov chain is reversible with respect to the distribution of G(n,p). We 
are now in a position to compute the terms needed to apply Theorems 13.61 
andO 

Lemma 4.15. Let (W, W) be the number of triangles in the exchangeable 
pair of Erdos-Renyi graphs (G,G') as defined above. If Q\(G) = W[W' = 
W + 1\G], then 

VarQi(G) 

< ^7^/(i - p)(i - ? 2 r~ 3 (i - p 2 (i - p)(i - p 2 r~ 3 ) (4.i7) 

A( n ~ 2 ) 

+ -W^a "P) 2 ((l " V +P 3 ) n - 4 -p(l -P 2 f n - G ) (4.18) 
, ^2 pf(i-pifn-*(l -p- p (l- p 2 ) 2 ) (4.19) 



p 6(l _ p ) 2 ((l _ p)n-3 (1 + p _ p2 )n -5 _ (1 _ p2) 



12(«- 2 ) 



2\2n-6 s < 



(5) 

(4.20) 

+ ^CAp 6 (1 " P) 2 (l " P 2 ) 2n " 9 (-2p + 4p 2 - 3p 4 + p 6 ) (4.21) 
I2J 
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+ A^-V(l -P?(l -p 2 ) 2n " 10 (V - 7p 4 + 4p 6 - p 8 ) (4.22) 
+ ^£p-p 6 (l -P) 2 (! -p 2 ) 2n - 10 (V - 7p 4 + 4p 6 -p 8 ). (4.23) 



Proof. Let JQj be the indicator that there is an edge between vertices i and 
X 

{i,j,k} with elbow i. We easily find 



j and V^ CJ := XijX^f. be the indicator that there is a T^-star on the vertices 



where we define the indicator variables 

From this point we note that the variance of Q\[G) is a sum of covariance 
terms times p 2 /i^) ■ For fixed and k, there are 3 Q) terms of the form 
Cov(YP' k ,Yu' ), where we are including Vax(Y?' ). In order to compute this 
sum, we will group these covariance terms with respect to the number of 
indices Y?' k and Yu '* share, which will yield the lemma after computing their 
covariances. 

As an example of the type of calculation involved in computing these 
covariance terms, note that 

E^' fc =p 2 (l-p)(l-p 2 ) 2n - 6 , 

so that 

Var Y?' k = p 2 {l -p){l- p 2 ) 2n " 6 (l - p 2 (l - p)(l - p 2 ) 2n " 6 ). 

Furthermore, for j ^ s, we find that 

E{ Y >> k Y?> k } = p 3 (l - p) 2 ((l - pf + 3p(l - p) 2 + P 2 (1- p)Y~\ 

Below we focus on carefully spelling out the number and types of covariance 
terms that contribute to the variance of Q\{G) and leave to the reader 
detailed calculations similar to those above. 

If {i,j,k} = {u,8,t} and Y/' k ^ Y°'\ then E{lf fc Y u M } = 0, so the 
corresponding covariance term is negative, which we bound above by zero. 
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In the case that if* = Y°> 1 , we obtain a variance term which corresponds 
to ()4.17p in our bound. 

Assume now that k} and {u, s, t} have exactly two elements in com- 
mon and consider which indices are equal. In the cases that Y^' is equal 
to Yi' k , Y?'*, Yu' k , Yu' j or Y^\ then Efl?'**^'*} = 0, so the correspond- 
ing covariance term is negative, which we bound above by zero. The two 
remaining cases to consider are Yu^ = Y^ ,k which contribute 2(n — 3) equal 
covariance terms leading to (|4.18p . and Yu = Yj' which also contribute 
2(n — 3) equal covariance terms leading to (|4. 19j) . 

Assume {i,j,k} and {u,s,t} have exactly one element in common; we 
have four cases to consider. There are 2(n — 3)(n — 4) covariance terms of 
the basic form Yu' 1 = Yi' 1 which leads to (I4.20p . there are 2(n — 3)(n — 
4) covariance terms of the basic forms Yu = Yu or Yu' 1 = which 

leads to (14.2ip . and there are f" - ^ 3 ) t erms of the form Yu 1 = Y^ 1 which 
yields (j02|) . 

Finally, if {i,j, k} and {u, s, t] are distinct sets, of which we have 3( n 3 3 ) 
ways of obtaining l^f'*, the corresponding covariance terms contribute (|4.23p 
to the bound. □ 

Lemma 4.16. Let (W, W) be the number of triangles in the exchangeable 
pair of Erdds-Renyi graphs (G,G') as defined above. If Q^\(G) = 1P[W' = 
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W -l\G], then 

Var Q_!(G) < ^t=A? 3 (1 - P?(l - p 2 ) n " 3 (l - p\l - p 2 ) n ~ 3 ) 

+ ^y^p 3 (i -P) 2 ((i - V +P 3 )- 3 -P 3 (i -p 2 ) 2 "" 6 ) 

+ ^&-W - p? ((1 - p)(i - 2 P + p 3 )"- 4 - P (i - p 2 ) 2 "" 6 ) 
+ - p) 2 (i - p 2 ) 2n ~ 8 (1 - p - Ki - p 2 ) 2 ) 

+ — kfVa - p) 2 ((i - p) 3 (i - 2p + p 3 ) n " 5 - (1 - p 2 ) 2n " 6 ) 
+ ^G-^ 1 - p) 2 u - p 2 ) 2n - 9 ((i - p) 2 - (i - v 2 ?) 

+ ^tP-/(1 " P) 2 (l " P 2 ) 2n ~ 10 ( 4p 3 - 7p 4 + 4p 6 - p 8 ) 

+ ^GfVu - P) 2 U - P 2 ) 2n " 10 (4p 3 - 7p 4 + 4p 6 - p 8 ). 

Proof. As in the proof of Lemma l4.15t let j be the indicator that there is 
an edge between vertices i and j and := X^jX^ be the indicator that 
there is a V-star on the vertices k} with elbow i. We easily find 

q-i(g) = ^ E E IK 1 - 

12; {j,k}i^j,k i&,j,k 

and from this point, the proof is very similar to the proof of Lemma f4.15l □ 

We are now in the position to prove Lemma 14.141 

Proof of Lemma \4-14\ Recall the notation of Lemmas 14.151 and 14.161 We 
have the following easy facts: 

(i) ~EW = Qp 3 , 

(it) a\ = Var^= rj(p 3 (l-p 3 ) + 3(n-3)p 5 (l-p)), 
(Hi) q 1 = MQ 1 (G) = (n-2)p 3 (l-p)(l-p 2 )"- 3 . 
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The first bound in (|4.16p now follows from Theorem 13.61 after noting that 
for (W, W') the number of triangles in the exchangeable pair of Erdos-Renyi 
graphs (G, G') as defined above and n a p — > c > for 1/2 ^ a < 1, then 

Var Qi{G) = O (p 5 ) , Var Q-\(G) = O (p 3 /n) , 

which follows easily from Lemmas 14.151 and 14.161 above. 

In order to prove the second bound in (|4,16p we will apply Theorem 13.71 
with G = G(n,p) an Erdos-Renyi random graph, G' obtained by taking 
a step from G in the Markov chain (reversible with respect to the law of 
G(n,p)) defined previously, and G" obtained as a step from G' in the same 
Markov chain. Setting (W, W, W") to be the number of triangles in the 
graphs (G,G' ,G"), and defining (as per Theorem 13. 7p 

Qi,i(G) = P[W" = W + i, W = W + i\G], 

it is easy to see that 

\2J {j,k}ij£j,k {s,t}^{j,k}u^s,t 

where as in the proof of Lemma 14.151 we define 

Y?' = (1 - X j)k )XijX i)k [ (1 - XijXik). 

From this point we find 

M\Q hl {G)-Q l (Gf\ = ^- qu 

since for fixed {j, k}, only one of the set {Y^ ,k }^ =1 can be non-zero. A similar 
analysis shows 

E\Q^ 1 (G)-Q^ 1 (G) 2 \ = ^-^ qi , 

and the second bound in (|4.16j) now follows from Theorem 13. 71 after collecting 
the pertinent facts above. □ 
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